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Abstract 

Perfect reconstruction filter banks can be used to generate a variety of wavelet bases. Using IIR 
linear phase filters one can obtain symmetry properties for the wavelet and scaling functions. 

In this paper we describe all possible IIR linear phase filters generating symmetric wavelets 
with any prescribed number of vanishing moments. In analogy with the well known FIR case, we 
construct and study a new family of wavelets obtained by considering maximal number of vanishing 
moments for each fixed order of the IIR filter. Explicit expressions for the coefficients of numerator, 
denominator, zeroes, and poles are presented. 

This new parameterization allows one to design linear phase quadrature mirror filters with many 
other properties of interest such as filters that have any preassigned set of zeroes in the stopband 
or that satisfy an almost interpolating property. 

Using Beylkin's approach, it is indicated how to implement these IIR filters not as recursive 
filters but as FIR filters. 
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1 Introduction 

Quadrature mirror filters (QMFsJl] allow one to design two-channel perfect reconstruction filter banks 
and therefore to generate wavelet bases. For implementation reasons we either consider finite impulse 
response (FIR) or realizable infinite impulse response (IIR) filters. The first class corresponds to 
polynomial filters, the second, wider class, to rational filters. In both cases, we consider filters with real 
coefficients and we identify filters with their z-transforms. Rational QMFs admit properties difficult 
or impossible to achieve with polynomial QMFs. For example, neither linear phase (except for the 
Haar filter) nor interpolation properties can be exactly obtained even though both properties can be 
approximately achieved using FIR coiflets |14j . In contrast, both properties can be exactly obtained 
with IIR filters but not simultaneously. Parameterizations of polynomials \15\ [22j or rational QMFs 
[9] are well known. In both cases, linear conditions on the coefficients of the filter, such as vanishing 
moments for the wavelet or the scaling function, become non-linear conditions on the parameters 
describing the coefficients. This problem can sometimes be overcome, and there are explicit families 
of rational filters with an arbitrary number of zeroes at —1, that is an arbitrary number of vanishing 
moments for the wavelet. But the additional requirement of symmetry conditions, which correspond 
to linear phase filters, again forces us to solve non-linear systems. For this reason, there are only some 
examples of linear phase filters generating wavelets with vanishing moments [10^ [T8] . 

Often, it is much simpler to describe the filters in terms of their absolute values (as in [3, 6j for 
polynomials or ^lOj for rational QMFs), leading to implicit constructions. In most cases the solutions 
are obtained numerically and the exact coefficients are not known except for small filter lengths. 
These numerical solutions are QMFs within a certain accuracy. In practical applications this is not 
an inconvenience because we always work within a certain precision that can in principle be attained 
except for large filter lengths. Nevertheless, we may like to incorporate the accuracy sought as a 
parameter in the process of design-implementation of the QMF. 

Leaving aside the problem of accuracy once the size of the filter increases, we come across a much 

deeper problem. Namely, how does one simultaneously impose different properties in an exact or 

approximate sense? A good example is Daubechies' parameterization of the magnitude response of 

polynomial QMFs. Even though she used it to obtain the autocorrelation coefficients of maximally flat 

filters, it is not clear how to use this description to obtain other designs. For example, when obtaining 

coiflets, neither her method [8j nor the simpler description in [14] can prevent bad behavior of the 

^Originally, filters leading to perfect reconstruction were named conjugate quadrature filters [20] but we use the term 
QMF as in 21 Pages 162-163]. 
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frequency response or guarantee zeroes in locations different than vr. Theoretical results are also hard 
to obtain, and we cannot even establish the existence of coiflets for an arbitrary number of vanishing 
moments. 

Here we propose a different approach for the design and implementation of QMFs. We first isolate 
a specific property of interest, like interpolation or linear phase, and then parameterize all possible 
rational QMFs with that property. The parameters can be used to impose additional properties or to 
optimize a particular design (13] . 

In this paper we choose linear phase as the initial property in our design. It implies two possible 
symmetries for the wavelet and scaling functions. Either the scaling function is even (and the wavelet 
is even about 1/2) or it is even about 1/2 (and the wavelet is odd about 1/2). These two situations 
give rise to two different classes of linear phase filters that we denote 0-SYM and 1/2-SYM. In this 
paper we discuss the 0-SYM case. For the 1/2-SYM case see [TOl ISl [13] . Even though there are no 
linear phase interpolating filters, it is possible to obtain linear phase and almost interpolating filters. 
It is enough to ask for a symmetric wavelet with enough vanishing moments. Both properties imply 
that the scaling function has vanishing moments (twice the number for the wavelet) and that assures 
the almost interpolating property }14] . 

Since linear phase filters have poles inside and outside the unit circle, we indicate how to implement 
these IIR by approximating them by cascades of FIR filters as in [1]. The FIR approximation does 
satisfy the quadrature mirror condition with any desired accuracy. The factors in the approximation 
directly depend on the location of the poles of the IIR filter, what allow one to introduce the accuracy 
as part of the filter design process. 

The main result of this paper is the complete description, with closed form expressions, of all 
possible 0-SYM filters generating wavelets and scaling functions with symmetry properties and any 
prescribed number of vanishing moments. Each 0-SYM filter is described using a minimal set of 
parameters which are directly related with the zeroes of the filter. This property allow us to construct 
filters with any preassigned set of zeroes in the stopband. In particular, we construct, for each fixed 
order, the filters which generate wavelets with maximal number of vanishing moments. We also give 
sufficient conditions on the parameters to assure that the frequency responses of the filters are always 
positive or that the filters have all their poles on the imaginary axis. 

The paper is organized as follows. In the next section we state some basic results and definitions 
and then identify the two kinds of symmetry. In Section [3] we describe, in closed form, all solutions 
of the 0-SYM case. We then identify a minimal set of parameters defining each solution, use them 
to obtain a variety of properties, and show explicit examples. The properties of the maximally fiat 
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0-SYM filters and an example of its FIR implementation are discussed in Section [H 



2 Preliminaries 

• Unless otherwise indicated, x and ^ are real variables while z is a complex variable. 

• A QMF is a 27r-periodic function mo, 

"^o(e) = E^'^«"'^' (2.1) 

kez 

such that 

|mo(e)P + |mo(e + ^)|' = l. (2.2) 

The numbers {/ifc} are the coefficients of the filter niQ. We assume rriQ to be a rational real 
trigonometric function. That is, {/ifc} are real and have exponential decay, but maybe an infinite 
number of them is non-zero. Condition (|2.2p is then equivalent to 

hkhk+2n = Sno for n G Z, (2.3) 

k 

where 6nm is defined as 6nm = 1 if m = n, and 6nm = otherwise. 

We denote by H the z-transform of {/i^}, H{z) = h^z^. We also refer to such H as a QMF. 
Its frequency response is H{e~^^) = mo(0- 

From (|2.3p . satisfies the functional equation (the QMF equation), 

H{z)H{z-^) + H{-z)H{-z-'^) = l. (2.4) 

In order to generate a regular Multiresolution Analysis (see [1]), we need two additional proper- 
ties for H. 

The first one is the normalization or low-pass condition. It forces 

mo(0) = 1 for H{1) = 1. (2.5) 
The second one assures that H is non-zero in certain locations of the unit circle [5]. 

Cohen's condition There exist no non-trivial invariant cycles for ^ — )• 2^ such that |mo| = 1 
on the whole cycle, i. e., there is no a; 7^ with 2"a; = oo (mod 27r) such that for all < j < n 

mo(2^a; + tt) = 0. (2.6) 

In practice, we first find a normalized H satisfying the QMF equation, and then verify Cohen's 
condition. 
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• A solution (p of 

v{^) = 2Y,hw{x-k) (2.7) 

k 

is called a "scaling function". Equivalently, 

(^(20 = mo(e)^(e), (2.8) 

where (p{i) = ip{x)e~^^^dx and mo is the associated filter. 
From (1231) and UO^ . 

CO 

^iO = n "^0(2-^^0 and 0(0) = 1. (2.9) 
k=l 

The mother wavelet ip is defined on the Fourier side as 

4^{2^) = mi{Om (2.10) 

where mi is a 27r periodic function, 

mi(e) = e-'«A(e)mo(C + vr). (2.11) 
Here A is vr-periodic and |A(^)| = 1 a.e. and mo is the complex conjugate of mo. 

2.1 The moment condition 

One fundamental property for wavelet bases [SJ \T2\ [T3] is the property of vanishing moments of the 
wavelet or the scaling function: 

/ x''i;{x) dx = for < k < M, (2.12) 
Jr 

/ x''ip{x) dx = dko for 0<k < N. (2.13) 

JR 

In terms of the filter they become [3 Page 258] 

D^H{-1) = for < A; < M, (2.14) 

D^H{1) = 6ko for < /c < iV. (2.15) 

2.2 The symmetry condition 

We now find which filters generate symmetry conditions on the scaling and wavelet functions. For the 
scaling function to be symmetric, we need real constants A and c such that for all x in R ip{X — x) = 
cip{X + x), that is 

'f{-0 = ce^'^^m- (2.16) 



By (|2.9p c should be 1 while using (|'2.8|) and (j2.16p with ^ and 2^, it follows that symmetry for 93 is 
equivalent to 

mo{-0 = e2'^«mo(e). (2.17) 



Since tuq has real coefficients, ?tio(— = rriQ^^) and we can write mo(0 = a(^)e '"^^j where 
a(^) = e^'^^mo(^) takes real values. In conclusion, symmetry for 99 is equivalent with linear phase for 

UlQ. 

In terms of H, (j2.17p becomes 

= z'^^H{z~'^). (2.18) 

Since replacing ip by an integer translate still generates the same Multiresolution Analysis, we can 
consider A in the interval [0, 1), and then the only possible values for A are or 1/2. 

Thus, either the scaling function is even (0-symmetry) or it is even about 1/2 (1/2-symmetry). 
By (12. 4p and (I2.18P the corresponding conditions on the QMF are 

0-SYM 

H^{z)+H^{-z) = l (2.19) 

H{z) = H{z-^). (2.20) 

1/2-SYM 

H^{z) - H^{-z) = z (2.21) 

H{z) = zH{z-^). (2.22) 

We now show which kind of symmetry we impose on the wavelet if we start with a scaling function 
with symmetry \ = Sip. 
As in ()2.16p . we require 

^(-0 = ce2-v.«^(e), (2.23) 

and then c is ±1 because '0 has norm equal to 1. 
Using (j2lll and (f2lT]) . 

i,{2i) = mo(e + 7r)e-^«d(20^(0 

where d(^) is a rational trigonometric function and (i(^)d(— ^) = 1. 
Because of (l2T6]) and (12171) 



7, _ e^gmo(-g-7r)(^(-g) _ e^'ge-^"'^^V^(20 
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Hence, to obtain (j2.23|) we need d{^) = zte'™"^ with m E Z and 2s^ = 1 — 2m . Choosing d{^) = — 1, 
that is 



"^i(e) = -e-'«mo(e + 7r), (2.24) 

we obtain = 1/2 , unique up to integer translations. Consequently, V is even about 1/2 or odd 
about 1/2 depending on whether s^p is or 1/2. An example of the latter case is the Haar basis, where 
the wavelet and scaling function are = Ijq 2^ — 1[2 and ip = l[o,i)- Here li is the indicator function 
of the interval I. 

The choice of mi in (I2.24p can be written in terms of the z-transforms of the filters as 

G{z) = -zH{-z-^). (2.25) 

3 The characterization of the 0-SYM case 

First, we introduce some additional notations and definitions. 

Sometimes we omit the variables of a function as in the following transformations 

/ = f{z) = f{-z) (3.1) 

r = nz) = f{z-'). (3.2) 

We use the notation s{ for the sign of the function / at the value t, 

4 = sign(f(t)) (3.3) 

If is a polynomial of degree d its reciprocal is z'^E{z~^). 
We call a an all-pass function if 

a{z)a{z-^) = 1. (3.4) 
Rational all-pass functions can be described as: 

a{z) = (3.5) 

where n is an integer, is a polynomial coprime with its reciprocal and E does not vanish at zero. 

The order of a rational function ^, where P and Q are coprime polynomials, is the maximum of 
their degrees. 

We use o for composition of functions. If n is a positive integer, /'"l denotes / o . . . o / (n times), 
/M is the identity and /t""*^^ is the inverse of /. 
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The bilinear transformation 



m = (3.6) 

1 + z 



maps the unit disk onto the right half plane {z : > 0} and the unit circle onto the imaginary 



/3(e^«) = -itan(|). (3.7) 



axis: 



The transformation /? gives a natural way to relate analog with digital filter design; see [21j and 
[161 Page 219], for more details. The conjugation by /3 of / is //3 = /3 o / o /j'"^! and similarly for 
/3(z)=/3(-z). 

The next theorem provides an explicit description of all 0-SYM filters. Conjugation by (3 is the 
key element for the simplicity of this parameterization. In Corollary [2] below we explain how to obtain 
vanishing moments for the associated wavelet and scaling functions. 

Theorem 1 There is a bijection between the set of 0-SYM filters H and the set of even all-pass 
functions a given by: 

V2H{-z) ' 

IJ H{\) — 1, H can be described as 

Hiz) = Mz){A{z)±V2Ai-z)) 

^ ' A{zY + A{-zY ±^/2A{z)A{-zy ^ ' ' 

where A{z) is a polynomial of even degree, coprime with A{—z), and equal to its reciprocal. Also, 
A{0)A{1) / and A{-1) = 0. 

Proof The transformation /3 satisfies: 

/3[-i]=/3, (3.11) 



Conjugation by /? yields, 



r = -/3, (3.12) 
/3 = ^. (3.13) 



d is all-pass <;=^ is an odd function, (3.14) 
d = d* <^=^ di3 is an even function. (3.15) 
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viz) = no — = — . (3.16) 



Define the function v as 



First we show that the maps in (j3.8p and (j3.9p are well defined. We need to check that for a even 
all-pass 

H = uoa{z'^)o(3 = uol3o(-)0 (3.17) 

a 

belongs to 0-SYM. Because a is even, o/j equals {ap) *, which gives ()2.20p . On the other hand, (I2.19P 
is invariant under odd functions and is odd, hence we only need to verify that v o j3 satisfies 
(|2.19p . This follows replacing z by [i{z) in the identity 

v'^{z) + v'^{z-^) = 1. (3.18) 

Next we start with H with 0-SYM and check that a in (|3.9p is even and all-pass. Using (|3.14p 
and ()3.15p it suffices to show that 

l + H-H 
P ° 1=-= 

is odd and invariant under z^^. The latter follows from (I2.20p . while using (I3.12|) the former is 
equivalent to the identity 

l + H-H V2H 
p o — = p o 



or 

\2 



V2H l + H-H 

1-{H -Hf = 2HH 



which follows from (|2.19p . 

To establish that the maps in (j3.8p and p.9p are each others inverses, let us denote them by S 
and T respectively. 

Let a be an even all-pass function. Since 

5(a) = u o a* o p 

we have 

l+S{a) -S{a){-z) _ 1 + t/(a*) - t/(a) 
V2Sia)i-z) ~^ V2i^{a) ^° 
Consequently, using that a is all-pass and the identity 

^ = 

we have 

ToS(a) = {^-^ — -){a) = a. 
9 



Finally, if if is a 0-SYM filter, 



l + ^/2r(i/) = (^)o/3, 



SoT{H) = poC- — ^ = _ — ^ L- 

H^ + {l-Hf 
I- H'^ - {1 + H"^ -2H) _ ^ 
1-2H + {H'^ + iJ2) 
For the second part of the theorem, we write the even all-pass function a as 

a(.)=.,n-ir+^-,^jgL_ (3.19) 

where r > 0, m > 1 (because H{1) = 1), B{z) = X^j=o^j-^"' '^ith bobr / 0, and B coprime with its 
reciprocal. In particular B{1) ^ 0. 
Thus, 

a*{f^{z)) = s-{-ir+^^t^ (3.20) 

where 

A{z) = (l + z)2™+2'^S(/52(z)) (3.21) 

r 

= (l + z)2"^6j(l-z)2^(l + z)2'-2i. (3.22) 

j=0 

The properties for A follow from the conditions on B while ()3.10p is the result of applying v to (j3.20p . 
■ 

As stated in the introduction, we can obtain an almost interpolating 0-SYM filter by simply 
choosing an associated all-pass function with a high number of poles at 0. 

Corollary 2 Consider a and H as in the previous theorem and let m > 0. The following conditions 
are equivalent: 

1 The multiplicity of as a pole of a is 2m. 

2 The multiplicity of —1 as a zero of H is 2m. 

3 The multiplicity of 1 as a zero of \ — H is 4m. 

// any of these conditions holds then by ()2.14p and (j2.15p 

x'^^(x) dx = Q for < A; < 2m, and 

I x^ip{x) dx = 6j;Q for < k < 4m, 
Jr. 

where ip and ip are the wavelet and scaling functions associated to H . 

10 



R 



Proof We simply show that the three conditions are equivalent with the polynomial A of (|3.22|) 
having 2m zeroes at —1. Part 1 follows from ()3.19p and ()3.22p while (j3.10p gives Part 2 (note that 
A{1) ^ but A{1) = because m > 0). The equivalence with Part 3 holds because (jS.lOp yields 

1 _ f{(2:) = M-z) /o 

A{z)'^ + A{-z)^±V2A{z)A{-z) ^' ' 

where the polynomials in the numerator and denominator are coprime. ■ 

Remarks 3 1. We can use the theorem as a dictionary to reformulate properties of H in terms of 
a and vice versa. Let us list some of them: 

zoG^-Hl} ^ a(/3(zo)) = oo (3.24) 
zoG^-^-l} ^ a(/3(zo)) = --^ (3.25) 
zoGF-1{oo} ^ a(/3(zo)) = e^^"^/^ (3.26) 

Because of symmetry and the QMF condition, the frequency response of H only takes values 
in the real interval [—1, 1]. In Table [1] we listed the correspondence between values of H on the 
unit circle and values of a on the imaginary axis. This correspondence can be easily obtained 
from ([3T711 . 

2. Because of (|2.19p . the zeroes of H are the negatives of the preimages hy H oi 1 and —1. These 
two sets of preimages are nicely represented in the numerator of (j3.10p . The zeroes of A{—z) 
(counted twice) are exactly the set H~'^{1} (see (|3.23p ) while the set H^^{—1} is given by the 
zeroes of A{—z) ± ^/2A{z) (counted twice) because 

iAi-z) ± V2A{z)f 
l + H(z) = ^ / ^' . 3.27 

^ ' A{z)^ + A{-z)^±V2Aiz)A{-z) ^ ' 

The even multiplicity of the preimages of 1 and —1 agrees with the fact that the frequency 
response of H cannot take values outside [—1,1]. 

3. The filter H cannot have poles either on the unit circle or on the real line, because /3 maps those 
values to the imaginary and real axes, and o, being even and real, maps them to the real line. 

Also, if H were a polynomial, H^^{(X)^ = {oo}, and we would have 



a(-l) 



,±i7r5/4 



leading to a contradiction because a has real coefficients. Therefore, we recover the result on 
the absence of real polynomial QMFs generating wavelets with even symmetry [20^ [6] . 
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m+r a 



4. Using the notation in (|3.3p . it follows from (|3.19p that = (—1) 
Since H{0) and H{i) are real numbers, (|2.19p implies 

H H 



Finally, ([3:9]) gives 



H{0) = and H{i) = mo(7r/2) = (3.28) 

\/2 v 2 



= and s'} = sf. (3.29) 



H{e'^) a(-itan(|)) 



(-1,0) (-V2,-;^)U(-^,0) 
{0} {-^/2,0} 
(0,1) (-oo,-V2)U(0,+oo) 
{1} {00} 

Table 1: Correspondence between values of H and a 
We now identify a minimal set of parameters defining each 0-SYM filter. 

Parameterization in terms of preimages of one 

Assume H{1) = 1. From (I3.24p . the set of preimages of 1 counted with their multiplicity (the set 
H^^{1}) completely determines the poles of the all-pass function a. Therefore, up to a sign, a, and 
thus H, are completely determined by H~^{1}. We fix the sign by choosing sf, that is sf^ . 

The conclusion is that two solutions of 0-SYM are the same if their sets of preimages of 1 and their 
values at i are identical. 

We now study the degrees of freedom in the set H^^{1}. 

First, because of (I3.23p . H^^{1} is the set of zeroes of A{—z)'^. Using the properties of A in (I3.10p 
and (|3.22p . A{—z) has 2(m + r) zeroes, 2m of them being 1. Also, if w belongs to H~^{1}, w and w^"^ 
belong to H^^{1} but —w does not. 

Therefore, once and m are fixed, there are only r degrees of freedom for a 0-SYM filter of order 
4(m + r). We isolate these r parameters in a subset A of H~^{1} by discarding reciprocals and the 
value 1, 

A = {Ai,---,AJ. (3.30) 
12 



In the next theorem, we rewrite our previous parameterization of H to show the dependence on A. 



Theorem 4 Let H be any 0-SYM filter with 2m zeroes at —1 and A the subset of H ^{1} defined 
in (j3.30|) . If is the function defined in (j3.16p and rj is the function 

riiz) = (3.31) 

then 

H{z) = V l^.f /3(r?(z)r n/3(-^)j . (3.32) 

Proof Since (ICTTl) and (IHr^Uj) imply 

H{z) = ^(a(i)) = ^.(.^(-1)-+^'^M__!) ), 

it suffices to show that 

s\{-\r^^^^ = u(r^{£)) (3.33) 

where 

r 

u(z) = .^/3(z)-n/5(^)- (3-34) 
J=i ^(^i) 

We know that A has degree 2(m + r), it is equal to its reciprocal, and (from ()3.23p ) its zeroes are 
the negatives of the elements of the set A considered together with their inverses. Therefore, there 
exist a constant c and a polynomial i?, 

r 

i^(z) = [](z + r?(A,)), 

such that 

A{z) = cz'"+'^(l + r,{z)rE{r]{z)) 

and the theorem follows. ■ 
We have now a straightforward way of designing a 0-SYM filter H by imposing the locations where 

the filter takes the value one. 

Since the value mQ(^) gives an idea of how narrow the transition band is, we will express it in 

terms of H~^{1}. The next corollary implies that if we trade vanishing moments for any other value 

on the arc (§,7r) of the unit circle (the stopband), then the slope of mo is steeper at ^. 



Corollary 5 Let H as in Theorem^ and ?7io(^) = H{e *^). Th 



en 



m'o(|) = -(2 - V2sf){m + (3.35) 

3 



2' - ^^(A,.). 
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Proof From the proof of Theorem HI we have 

moiO = i^(n(cos(C))). (3.36) 
Taking derivatives at vr and using n(0) = we obtain 

) = -.'(sf K(0) = -^^^^n'(0). 

Since 

m 

the logarithmic derivative of ()3.34p at is, 

The next two results are obtained imposing conditions on the set H^^{1}. 

Proposition 6 Let H be a 0-SYM filter with = 1 and 
H-^{1} n {z : \z\ = 1} = {1}. Then 

H{e'^)>0 for Ce(-^,7r). 
Proof Since H{1) = 1 and H{z^^) = H{z) it suffices to consider ^ in (0, vr). From (I3.19p . 

aim)) = a{-^tM|)) = ^n^frsl^iy (3-37) 

where x = tan^(|), m > 1, and E{x) = B{—x). 

The zeroes of E are — /3(— A)^ where A 7^ 1 belongs to H^^{\}. (See (j3.2ip and the second remark 
in Remarks [31) These zeroes cannot be positive numbers because 

-/3(-A)2 G (0,+oo) ^ /3(-A) G Ri ^ |A| = 1, 

which contradicts our assumption on H~^{V\. 

Thus, E has constant sign in (0, +cxd) and (j3.37p implies that a[f3{e^)) is positive for ^ G (0,vr). 
The result follows using Table [TJ ■ 

We now discuss a typical filter design problem. Suppose the desired frequency response vanishes at 
certain specific locations of the stopband. By choosing appropriate preimages of 1 in the passband, we 
can obtain the given zeroes but, in principle, we cannot avoid further zeroes in the stopband because 
they are related to preimages of —1. The next theorem provides sufficient conditions to prevent this 
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last situation from happening and, incidentally, forces all poles of the filter to be purely imaginary. 
By restricting H~^{1} to the right half plane we can parameterize 0-SYM filters with any preassigned 
set of zeroes in the stopband. We can then obtain a filter with minimal order or, considering more 
parameters, even impose additional properties. 

Theorem 7 Let H be any 0-SYM filter with H~^{1} C {z : > 0}. We have the following 

properties: 

A H-^oo} C m, 

B H-^{-l} C {z : < 0}. 

Proof The assumption on H and (j3.24p imply that all poles of the corresponding all-pass function a 
are in the unit disk. Therefore, a{z~^) = is a finite Blaschke product and then a maps the unit 
circle into itself and {z : \z\ > 1} into the unit disk and vice versa. 

To verify Part A consider p G H'^{oo}. By (j3.26p . a{l3{p)) and then /3(p) belong to the unit circle. 
Thus, p is purely imaginary by (13. 7p . 

For Part B start with p G H^^{—1}. Since (I3.25P implies that a(J3{p)) belongs to the unit disk, 
we have |/3(p)| > 1 and then < 0. ■ 

Corollary 8 For {6*1 • • • 6* J m (f , tt) let 

where the set {Ai, • • • , A^} is contained in {z : ^^{z) > and \ z\ ^ 1} and satisfies the conditions listed 
for (I330D . 

Then, H is a 0-SYM filter that vanishes in the stopband only at the 2{m + s) values 

e9i, e~ 

Proof Using Theorem U define H with 2m zeroes at —1 and A = {— e^i, • • • , — eS^, Ai, • • • , At}. 
Theorem [7] assures that there are no other zeroes in the stopband. ■ 

Note that "wrong" choices of H~^{1} can lead to very poor frequency responses. In Figure [1] we 
plotted the frequency responses of two order 12 filters with = 1 and two zeroes at —1 but different 
choices for the other zeroes. Their set of preimages, zeroes, and poles are listed in Table [2j In one case 
we have chosen the zeroes in order to violate Cohen's condition. The poor frequency response is not 
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a peculiarity of this choice. Perturbing the values of the preimages we still obtain a similar response. 
Note that the zeroes either belong to the unit circle or they are real numbers. For the other example, 
all preimages belong to the right half plane and therefore all zeroes related to preimages of —1 belong 
to the same half plane. In agreement with the previous theorem, all poles are purely imaginary. 



1 

0.5 



-3 
■0.5 

-1 











I 


, 14 









3 . 14 




-3. 14 







3 . 14 



Figure 1: Two 0-SYM filters of order 12. On the left we violated Cohen's condition. On the right we 
choose all preimages of 1 in the right half plane. 



4 Maximally Flat Filters 

As in Daubechies' construction for compactly supported wavelets [6], we use our representation of 
the 0-SYM class to construct, for each fixed order of the filter, symmetric wavelets with maximum 
number of vanishing moments. For the reasons explained before, they correspond to the choice a{z) = 
s^(— in (|3.17p . For simplicity of notation, we write sf^ = (— 1)*^, where (5 is a positive integer. 

4.1 Description and Properties 

For 5 and n positive integers, we define 

Eliz) = (4.1) 



that is 



4nV / , , , / /TT/i 2\2n' ' 



(1 + + (1 - + (-1)^+'^ V2(l 



The filters Ef^ have maximum number of zeroes at —1 (2n in fact) in the set of rational QMF 
with 0-SYM and order at most 4n. For their frequency response, we use p.7p and ()4.ip : 
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K(-l)'tan2"|) 



K(;^ + (-l)^cot2-|)2- 



It readily follows that, as n goes to oo, Ef^ is approaching the ideal low-pass filter l 



(4.3) 
(4.4) 



(4.5) 



1 if < lel < tt/2, 
if tt/2 < 1^1 < vr. 

When 5 = the convergence is uniform because v is a decreasing function in the interval (0, 1). To 
illustrate the convergence, a few filters are simultaneously plotted in Figure [2] (for 5 = 0) and Figure 
[3] (for 6 = 1). 

Because of (12.25p . the filters mi associated with the wavelet converge to the ideal high-pass filter. 
Compare this situation with Daubechies' family of maximally fiat filters that converge to the ideal 
filters only in absolute value jllj . 

1 



.707 - 



. 4 



0.2- 




-3.14 3.14 

Figure 2: Approximating the ideal low-pass filter by maximally fiat filters -E^^- The filters shown in 
the picture correspond to n = 2, 3, 8, and 20. 
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0.5- 



-0 . 4 



■0 .707 - 




3 . 14 



Figure 3: Approximating the ideal low-pass filter by maximally flat filters El^. The filters shown in 
the picture correspond to n = 2, 3, 8, and 20. 

Next, we explicitly find all zeroes and poles of Ef^. In agreement with the results of the previous 
section, the frequency response of E^^ is positive in {—it, vr) (Proposition [6]), all poles of £'4^ are purely 
imaginary (Theorem [71 Part A), and except for —1, all zeroes of E^^ belong to the right half plane 
(Theorem [71 Part B). The latter result also holds for Daubechies' family of maximally flat filters [19^ 
Theorem 2]. 

Zeroes of Ef^ 



We already know that —1 is a zero of multiplicity 2n. The other zeroes are the values of uj such 
that H{-uj) = -1. Since a{z) = {-if+'^z-'^"' , using ([3:25]) we need to solve 



V2 



It follows that the zeroes are 



/3(2-^e^(2i+n+5-i) 



(4.6) 



for —n < j <n. Note that they belong to the right half-plane. 

With respect to zeroes on the unit circle, we now show that there is a pair of complex conjugate 
zeroes if 5 = 1 and no zeroes if 5 = 0. 
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We check in (|4.6|) when the argument of f3 is purely imaginary. Since —n < j < n, it foUows that 
sin(^(2j + (5 — 1)) = only if 6 = 1 and either j = oi j = n. In either case, the real part equals 

2(21^) = /3(2-2k). 
V 

For all n the numbers above are less than ^ and Cohen's condition holds. 



Note that for any real numbers A and 4) 

A ^ 77(A) + cos ^ 



pC4) = ^-^^^ (4.7) 



where rj is the function in (j3.3ip and 7(2) = — . 

Thus, except for —1, the zeroes of E^^ can be also expressed as 

/ 7(2^)-icos(^(2j+<^-l)) | 

for —n < j <n. 
Poles of 

By (|3.26|) , we need the values of z such that 

Because of (13. 7p and (13. lip the poles are 

(zbitan f — (5 + 8j + 46 + 4n) 
L V Ion 

for —n < j < n. 

Remark 9 The recurrence relation 

follows substituting z by ^ in ()4.ip and using the identity 



Thus, for all s, 



Ei^,,^ = Eir{{ihD, 

and the behavior of E^^ for large n is dictated by the iteration of the map 1/r/. This map is conjugate 



/3oio^[-il =z2. (4.9) 

7? 



to implying that its Julia set is the unit circle. We have then another explanation why these families 
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are approaching the ideal filter: When n goes to oo, restricted to the unit circle has only 

three possible values. In fact, 

' 1 if lei < ^/2, 

(l/77)[-](e^«) = <! oo if |e|=vr/2, (4.10) 

if 7r/2 < ICI < vr. 



Theorem [T] implies that 



Hg = /3oi?o/3[-^] 



1 



(1 + V2a)2 

for any 0-SYM filter H. When H{1) = 1, this conjugation changes the location of the fixed point from 
1 to 0. We can use /3 = (3{—z) to change it from 1 to oo 

H^{z) = {1 + V2aiz~^)f. 

The choice a{z) = (^—ly+^z^'^'"' yields the maximal family E^^ whose set of preimages of 1 equals 
{1, • • • , 1}. Thus, the conjugation by f3 of each Ef^ gives a polynomial! It is 



4n 



(1 + (-l)^+"V2z2'^)2. 



(4.11) 



0.5 



■3.14 



3 .14 



1 

0.5 

-3 
■0.5 



14 



3 .14 



For (5 = 0, 1 



1 + (-1)^+1 V2tan6(e/2) 



1 + (-l)'5+i^/2 tan6(^/2) + tani2(^/2) 
have maximal number of zeros at vr inside the set of rational linear phase QMFs of order 12. 

Figure 4: Frequency responses of the maximal filters E^2 (left) and EI2 (right). 



4.2 FIR implementation 

Following p], any IIR filter H{z) 



can be approximated on the unit circle by a FIR filter F{z) 



F{z) = z'^'Piz) H F, 



k\z 



(4.12) 



fc=0 
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The polynomials F]t[z) only depend on the zeroes of Q{z) and on the accuracy sought in the approx- 
imation. 

F{z) is not an exact QMF anymore, but if for w on the unit circle 

\H{w)-F{w) 



\H{w) 



then, for all z 



F{z)F{z-^) + F{-z)F{-z~^) = l + R{2 



where the FIR filter R{z) satisfies < 3e for w on the unit circle. We assumed e < fdegQ' 

As an example, we approximate E^2i^^^ positive maximally flat 0-SYM filter with six zeroes at 



— 1. (See Figure m) To avoid negative powers of z in (|4.12p . we shift the phase of the filter by the 
positive integer and approximate z^ Ei2{z)- 

The coefficients of the factors of the approximating filter for e = 10~^ are listed in Table [3l The 
factor Fq equals 1 because the denominator of is an even function. The wavelet and scaling 
functions generated with the approximating filter are plotted in Figures [5] and [6l 




Figure 5: Wavelet function for FIR approximation of the maximally flat filter Ei2- 
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Parameters 


Zeroes 


Poles 




-1 


-0.84955807 - 0.74802903 I 




-1 


-0.84955807 + 0.74802903 I 




-0.74212 - 0.67026705 I 


-0.66304573 - 0.58380642 I 




-0.74212 + 0.67026705 I 


-0.66304573 + 0.58380642 I 


sf = m=l 


-0.30901699 - 0.95105652 I 


-0.40197132 I 




-0.30901699 + 0.95105652 I 


0.40197132 I 




0.17142917 


-2.4877396 I 


A = {eh\el'''} 


0.65396257 - 0.75652691 I 


2.4877396 I 




0.65396257 + 0.75652691 I 


0.66304573 - 0.58380642 I 




0.80901699 - 0.58778525 I 


0.66304573 + 0.58380642 I 




0.80901699 + 0.58778525 I 


0.84955807 - 0.74802903 I 




5.8333128 


0.84955807 + 0.74802903 I 




-1 


-0.083442717 I 




-1 


0.083442717 I 




-0.79015501 - 0.61290705 I 


-0.57528543 I 




-0.79015501 + 0.61290705 I 


0.57528543 I 


sf = m=l 


-0.56208338 - 0.82708057 I 


-0.73702991 I 




-0.56208338 + 0.82708057 I 


0.73702991 I 




0.03560146 - 0.65573566 I 


-1.356797 I 


— |g0.2l7rj g0.3l7rj| 


0.03560146 + 0.65573566 I 


1.356797 I 




0.036837087 


-1.7382676 I 




0.082552825 - 1.5205228 I 


1.7382676 I 




0.082552825 + 1.5205228 I 


-11.984269 I 




27.146554 


11.984269 I 



Table 2: Zeroes and poles of filters constructed via their preimages of one. 
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Factor 


Coefficients 


Factor 


Coefficients 


P 


-0.0001011263580012439 


F2 


1.348299677989997e-7 




0.0029296875 




0.007308809891655256 




0.01818488314800746 




0.162081739736554 




0.0537109375 




0.6612186310836452 




0.1156706046299813 




0.162081739736554 




0.193359375 




0.007308809891655256 




0.2324912771600248 




1.348299677989997e-7 




0.193359375 








0.1156706046299813 




0.0001276629992294306 




0.0537109375 




0.03971627520745388 




0.01818488314800746 




0.920312123586633 




0.0029296875 




0.03971627520745388 




-0.0001011263580012439 




0.0001276629992294306 


Fi 


-0.00268082617584078 


Fi 


1.925310635034675e-8 




0.6429247852752233 




0.001585818961857287 




-4.433610674839401 




0.996828323570073 




8.58673343148004 




0.001585818961857287 




-4.433610674839401 




1.925310635034675e-8 




0.6429247852752233 








-0.00268082617584078 


F5 


2.492072168636633e-6 








0.999995015855663 








2.492072168636633e-6 



Table 3: FIR approximation of a shifted maximal filter with six zeroes at — l.The coefficients hsted 
correspond to the factors in Eq. (I4.12j) . In this example = 0, n = 5, and Fq = 1. 
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5 Conclusion 



The results described in this paper allows one to explicitly construct all possible linear phase IIR 
quadrature mirror filters, when the slope of the phase is zero. Each filter is described using a minimal 
set of parameters which are directly related with the zeroes of the filter. This property permits to 
construct filters with any preassigned set of zeroes in the stopband. In particular, we construct, for 
each fixed order, the maximally flat fllters which generate wavelets with maximal number of vanishing 
moments. We also give sufficient conditions on the parameters to assure that the frequency responses 
of the filters are always positive or that the filters have all their poles on the imaginary axis. We 
indicate how to implement these IIR filters by approximating them by cascades of FIR filters. The 
FIR approximation does satisfy the quadrature mirror condition with any desired accuracy. The 
factors in the approximation directly depend on the location of the poles of the IIR filter, what allow 
one to introduce the accuracy as part of the filter design process. 
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